Kinetics of the superconducting charge qubit in the presence of a quasiparticle 



o 
o 

(N 

wo: 
< 



R. M. Lutchyn 1 , L. I. Glazman 1 , and A. I. Larkin 1,2 
1 W.I. Fine Theoretical Physics Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA 
2 L. D. Landau Institute for Theoretical Physics, Moscow, 117940, Russia 
(Dated: February 6, 2008) 

We investigate the energy and phase relaxation of a superconducting qubit caused by a single 
quasiparticle. In our model, the qubit is an isolated system consisting of a small island (Cooper-pair 
box) and a larger superconductor (reservoir) connected with each other by a tunable Josephson 
junction. If such a system contains an odd number of electrons, then even at lowest temperatures a 
single quasiparticle is present in the qubit. Tunneling of a quasiparticle between the reservoir and 
the Cooper-pair box results in the relaxation of the qubit. We derive master equations governing 
the evolution of the qubit coherences and populations. We find that the kinetics of the qubit 
can be characterized by two time scales - quasiparticle escape time from the reservoir to the box 
rr 1 and quasiparticle relaxation time r. The former is determined by the dimensionless normal- 
state conductance g T of the Josephson junction and one-electron level spacing 8 r in the reservoir 
(ri n ~ g T 5 r ), and the latter is due to the electron-phonon interaction. We find that phase coherence 
is damped on the time scale of rr 1 . The qubit energy relaxation depends on the ratio of the two 
characteristic times r and r ; ~ and also on the ratio of temperature T to the Josephson energy E ,. 

PACS numbers: 03.67.Lx, 03.65.Yz, 74.50.+r, 85.25.Cp 
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I. INTRODUCTION 

Recent experiments with superconducting charge 
qubits demonstrated coherent oscillations between two 
charge states of a superconducting island, a so-called 
Cooper-pair box (CPB), in a single-electron devicei&2£. 
A device with a large superconducting gap A > E c > 
Ej^>T can be controlled with the gate voltage and 
magnetic flux, and has only one discrete degree of free- 
dom: the number of Cooper pairs in the box. (Here E c 
is the charging energy of the island, and E r is an effec- 
tive Josephson energy of its junctions with the reservoir; 
E j can be tuned by the flux.) The practical implemen- 
tation of superconducting qubits requires long coherence 
timea^. Although the contribution of the quasiparticles 
to the decoherence of the existing qubits is not the lead- 
ing ona^ii, it limits qubit operations on the fundamental 
level. This motivates us to study the qubit dynamics in 
the presence of a quasiparticle. 

In this paper we consider energy and phase relaxation 
in a superconducting charge qubit due to the presence of 
a single quasiparticle. Passing of a quasiparticle through 
the Josephson junction leads to the escape of the qubit 
out of a two-level system Hilbert space, and thus deter- 
mines the decay rate of coherent oscillations. In a previ- 
ous paper—, we evaluated rates of the elementary acts in- 
volving quasiparticle tunneling. In particular, we showed 
that the rate of tunneling into the CPB Ti n is determined 
by the dimensionless (in units of e 2 /h) conductance g T 
of the junction between the CPB and the reservoir, and 
level spacing in the reservoir T ln ~ g T <5 r /47r. The rate of 
tunneling out of CPB is r ou t ~ g T 6f,/4TT with 5b being 
level spacing in the box. Taking into account the dif- 
ference in the volumes (V r ^S> Vt>), one can notice that 
Tout Tin- Therefore, for a sufficiently small box the 
quasiparticle dwelling time in the CPB is very short, and 



the qubit spends most of the time in the "good" part 
of the Hilbert space. Nevertheless, the evolution of the 
qubit will be affected by quick "detours" the qubit takes 
outside that part of the Hilbert space. We demonstrate 
that even a single detour destroys the coherence of the 
qubit. Combined with the phonon-induced relaxation of 
the non-equilibrium quasiparticle in the reservoir, the de- 
tours also lead to the relaxation of the populations of the 
qubit states. We derive and solve the master equations 
for the dynamics of the qubit, which describe its relax- 
ation caused by a quasiparticle. 

The paper is organized as follows. We begin in Sec. [H] 
with a qualitative discussion and brief overview of the 
main results. In Sec. lIIII we derive the microscopic master 
equations for the dynamics of the qubit without quasipar- 
ticle relaxation. We solve these equations and discuss the 
kinetics of the coherences and populations in the Sec. II VI 
andlvl In Sec. I VII we incorporate mechanisms of quasi- 
particle relaxation into the master equations and solve 
them for fast and slow quasiparticle relaxation limits. 
Finally, in Sec. IVIll we summarize our main results. 



II. QUALITATIVE CONSIDERATIONS AND 
MAIN RESULTS 



Dynamics of the superconducting charge qubit is con- 
ventionally described by an effective Hamiltoniar. 9 



H qh = E C {N - N g f + H, 



(1) 



where E c is charging energy of the box, N is the 
charge of the CPB in units of one-electron charge e, 
and N g is the dimensionless gate voltage. The Hamil- 
tonian of Cooper-pair tunneling H r is defined as H , = 

(\ N + 2 > ( N \ + H.c), where E, is the Josephson en- 
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ergy. In the case of a large superconducting gap 

A > E c > Ej > T, 

the existence of quasiparticles is usually neglected, and 
the dynamics of the system is described by Hamiltonian 
JTJ with only one discrete degree of freedom, the excess 
number of Cooper pairs in the box. The qubit can first be 
prepared in a state with an integer number of pairs in the 
CPB, e.g., in state \N) with even N, and then tuned to 
the operating point by adjusting the gate voltage to the 
value N g = iV+1. The charge degeneracy between states 
\N) and |JV + 2) at this point is lifted by the Josephson 
tunneling, and the states of the qubit are described by 
the symmetric and antisymmetric superposition of the 
charge states, R = |jV)+l f +2> and ^ 
with energies 



V2 



V2 



E c - E J2 and lo + = E c + EJ2, 



(2) 



respectively. Other charge states have much higher en- 
ergy, and effectively the Cooper-pair box reduces to a 
two-level system. Being coherently excited by such tun- 
ing, the qubit oscillates between the states |+) and |— ) 
with the frequency defined by the Josephson energy Ej. 

The presence of a quasiparticle with a continuum exci- 
tation spectrum provides a channel for relaxation of the 
qubit. If the state \N) is prepared in equilibrium condi- 
tions, then the quasiparticle resides in the reservoir part£ 
of the qubit. Upon tuning of the qubit from state \N) 
to the operating point, a charge degeneracy point for the 
system is passed at N g = N + 1/2, see Fig. 1. (Here- 
after we assume equal superconducting gap energies in 
the reservoir and Cooper-pair box.) However, if tuning 
is performed fast enough, the quasiparticle remains in the 
reservoir^. 

The coherent charge oscillations at the operating point 
of the qubit continue until the particle finds its way into 
the CPB. On average, this occurs on a time scale of the 
order of (g T 8 r / Aii)~ x . There are several assumptions that 
allow for this estimate^. First, the quasiparticle level 
spacings 5b and S r in the CPB and reservoir, respectively, 
must be small compared with the temperature T, which 
defines the initial width of the energy distribution of the 
quasiparticle. Second, the fluctuations of the potential 
between the grains must exceed 5 r , see, e.g., Ref. 
Third, we neglected the difference between A and E c 
when including in the estimate the density of states and 
tunneling matrix elements of a quasiparticle at energy 
~ E c above the gap in the CPB. Under these conditions, 
the average time it takes the quasiparticle to leave the 
reservoir and enter the CPB is of the order of the inverse 
level width of a state in the reservoir with respect to 
leaving it through the junction of conductance g T 



g T $r 

ill 



(3) 



Once the quasiparticle enters the CPB, the charging 
energy that the qubit has at operation point is trans- 
formed into the kinetic energy of the quasiparticle, see 




FIG. 1: (color online). Energy of the Cooper-pair box as a 
function of dimensionless gate charge N g in units of e (solid 
line). Near the degeneracy point (N g = 1) Josephson coupling 
mixes charge states and modifies the energy of the CPB. The 
dashed line corresponds to the charging energy of the CPB 
with an unpaired electron in the box. At N g — 0.5, the tun- 
neling rate Ti n lifts the degeneracy between the ground state 
of the CPB (solid line) and a state with a single quasiparticle 
in CPB (dashed line). We assume equal superconducting gap 
energies in the reservoir and CPB. 



Fig. 2] The quasiparticle may escape the CPB leaving 
the qubit in the excited or ground state, see Fig. 2. The 
rates of escape into these states are different due to the 
difference of the kinetic energies available to the quasi- 
particle upon the escape and due to the energy depen- 
dence of the superconducting density of states u(e). If 
the qubit ends up in the excited state upon the escape, 
then only energy e ~ T is available for the quasiparticle, 
and v{T) ~ (5,7 1 (A/T) 1 / 2 (we used here the condition 
A ^> T). The corresponding escape rate is 



r, 



47T 



(4) 



If the qubit arrives in the ground state, then energy ~ E 7 
is available to the quasiparticle, and its density of states 
in the final state is v{E } ) ~ 5~ 1 (A/E J 1 / 2 ; the rate of 
escape to this state is 



(•5) 



These two rates are much higher than because 5b 3> 
S r , so detours of the quasiparticle to the CPB are short 
compared to the time quasiparticle spends in the reser- 
voir. Nevertheless, the typical time the quasiparticle 
spends in the CPB is much greater than the oscillation 
period of the qubit. Indeed, the ratio 





r l-> 

L out 



(6) 



is small: 5b/ A ~ 10~ 4 — 1CP 3 for any reasonable size 
of the CPB (we used here the Ambegaokar-Baratoff re- 
lation between E J} g T and A). The times of return of 
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FIG. 2: (color online). Schematic picture of the transitions 
between the qubit states in the presence of a quasiparticle in 
the reservoir, e.g. |+, E p ) <-» |o, Ek) ■ Having kinetic energy 
~ Ej the quasiparticle can emit a phonon. The corresponding 
state of the system is | — , E p +Ej). 



the quasiparticle back to the reservoir are randomly dis- 
tributed. The probability of the quasiparticle returning 
to the reservoir during times that are short compared to 
the oscillation period 2ir/Ej is of the order T out /Ej and 
is small (here we do not distinguish between T\^ t and 
r^t). Therefore, a single detour of the quasiparticle into 
the CPB destroys coherent oscillations of the qubit with 
overwhelming probability. Taking into account the rela- 
tion Tin <C r out , we find that the dephasing rate for the 
qubit, induced by the quasiparticle, is limited by the rate 
of quasiparticle tunneling into the CPB 



i-2 



(7) 



with T in of Eq. ©. 

Unlike the phase, the energy stored in the degrees 
of freedom described by the qubit Hamiltonian Q is 
not dissipated at the short time scale given by Eq. 0. 
We start analyzing the time evolution of the qubit en- 
ergy by considering the limit of infinitely slow quasipar- 
ticle relaxation (the latter typically is determined by the 
electron-phonon interactioni2iI4). If initially the system 
was prepared in the |+) state, then upon a single cy- 
cle of quasiparticle tunneling, the qubit ends up in the 
ground state w ith a sm all probability defined by the ratio 

r out/ r ott ~ V t / e .j- In other words, the qubit energy 
will randomly change over time between two values w+ 
and uj-, see Eq. (J5J). The portion of the time, the qubit 
spends in the ground state and the quasiparticle is ex- 
cited to the energy E d above the gap, is small ~ y/T/E ,. 

The portion of time that the quasiparticle spends in an 
excited state in the reservoir becomes important when we 
account for the phonon-induced relaxation of the quasi- 
particle. Having energy E T , the quasiparticle may emit a 
phonon at some rate 1/r and relax to a low-energy state. 
The relaxation of the quasiparticle will prevent further 




reexcitation of the qubit into |+) state and result in qubit 
energy relaxation. To find the energy relaxation rate of 
the qubit, we multiply the portion of time the quasipar- 
ticle spends in the excited state by the relaxation rate 
1/r: 



(8) 



This estimate is applicable if r 3> l/rj n , and many cycles 
occur before the energy is dissipated into the phonon 
bath. 

In the opposite case of fast relaxation r <C l/r^, the 
quasiparticle loses its energy the first time it gets it from 
the degrees of freedom of the Hamiltonian QJ. Therefore, 
in this case the qubit energy relaxation on average occurs 
on the time scale 




(9) 



For aluminum, a typical superconductor used for 
charge qubits, the quasiparticle relaxation time r is in- 
deed determined by the inelastic electron-phonon scatter- 
m gi&14, and a t low energies (e < E ,) it can be estimated 
as t > 10 (is. For a small mesoscopic superconductor 
this time is longer than the typical values of 1/Fi n , and 
Eq. (JSJ gives an adequate estimate for the qubit energy 
relaxation rate. 

A comparison of the phase relaxation time Q with 
even the shortest of the two energy relaxation times ©, 
indicates that the coherence is destroyed much earlier 
than the populations of the qubit states approach equi- 
librium. Therefore one may consider the decay of qubit 
coherence separately from the process of equilibration, 
which involves the electron-phonon interaction in addi- 
tion to the quasiparticle tunneling. In the rest of the 
paper, we derive and solve the master equations, which 
yield results discussed qualitatively in this section. 



III. DERIVATION OF THE MASTER 
EQUATIONS WITHOUT QUASIPARTICLE 
RELAXATION. 



The Hamiltonian of the entire system consists of the 
qubit Hamiltonian -ff q b, BCS Hamiltonians for the su- 
perconducting box and reservoir H„„ a and HZ„„ , respec- 
tively, and quasiparticle tunneling Hamiltonian V: 



H = Hr, 



(10) 



where H = H b 



■HZ 



ifqb and perturbation Hamil- 



tonian V takes into account only tunneling of quasipar- 
ticles V = H T —Hj. The tunneling Hamiltonian H T is 
defined as 



=£(' 



„t 



+ /i.e.), 



(11) 



kpa 
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where t kp is the tunneling matrix element, Cfe )0 -, Cp !<T are 
the annihilation operators for an electron in the state k, a 
in the CPB and state p, a in the superconducting reser- 
voir, respectively; H } is of the second order in tunneling 
amplitude^ 

Hj = | AO (N\H T —±—H T \N+2) (N+2\+H.c. 

E — tlQ 

(12) 

The matrix element (N\ H T ^j^H T \N + 2) is propor- 
tional to effective Josephson energy E J: and H T is de- 
fined in Eq. Ijlll) . Without quasiparticles Hamiltonian 
Hq reduces to Eq. and qubit dynamics can be de- 
scribed using the states |+) and |— ). In the presence of a 
quasiparticle, qubit phase space should be extended. Rel- 
evant states now are \+,E p ), \ — ,E p ), and \o,E k )- The 
first two states |±, E p ) correspond to qubit being in the 
excited (ground) state |±) and a quasiparticle residing in 

the reservoir with energy E p = y 1 ^ + A 2 : 

\±,E P ) = \±)®\E P ). 

The third state |o, E k ) describes the qubit in the "odd" 
state with N+ 1 electrons in the box, i.e., the qubit 
escapes outside of its two-level Hilbert space. Here 
E k = V£fe + A 2 is the energy of the quasiparticle in the 
box. Perturbation Hamiltonian V causes transitions be- 
tween the states \±,E p ) and |o, Ek), see Fig. |3 Note 
that V does not induce the transitions between \ + ,E p ) 
and | — , E p ). 

The evolution of the full density matrix of the system 
is described by Heisenberg equation of motion (h = 1): 

p I (t) = -i[V I (t),p I (t)], (13) 

where subscript / stands for the interaction representa- 
tion, i.e., Vj(t) = e iH ty e -iH t_ The iterative so i u tion 

of Eq. 111. 'ill yields for the matrix elements of the density 
matrix 

(s\p I (t)\s')=-i{s\[V I ,p(0)]\s') (14) 

- fdrisUVMiVAt-rlpAt-rMs'}, 
Jo 

where \s) can be \ + ,E p ), \ — ,E p ), or \o,Ek). The interac- 
tion Hamiltonian V has no diagonal elements in the rep- 
resentation for which Hq and p(0) are diagonal. There- 
fore, the first term in right-hand side of Eq. I|14|) is equal 
to zero, 

(slpM^-fdTisliVMdVAt-T^pAt-TMs'), 

Jo 

(15) 

Equation (|T3|) implies that evolution of the projected 
density matrix is proportional to V 2 . Since interaction 
is assumed to be weak, the rate of change of p c (t — r) 
is slow compared to that of Vj(t). Therefore, one can 
approximate p r (t — r) by p I (t) in the right-hand side of 



Eq. (|15fl (see, for example, Refs. [l5llT^ | for more details 
on the derivation). Finally, going back to the original 
representation, we arrive at the following system of mas- 
ter equations 
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(16) 

where states \m), \n), \s), and |s') denote |+, E p ), | — , E p ), 
or |o, Ek), and the sum runs over all possible configura- 
tions. The system of equations (|16|) describes the kinet- 
ics of the qubit in the presence of a quasiparticle in the 
Markovian approximation. We are interested in elements 
of the density matrix that are diagonal in quasiparticle 
subspace, e.g., P + _(E p ,t) = (+,E p \ p(t)\E p ,—), since at 
the end one should take the trace over quasiparticle de- 
grees of freedom to obtain observable quantities. Note 
that Eq. i|16[l which describes evolution of a closed sys- 
tem (the qubit and the quasiparticle) does conserve its 
total energy We will include the mechanisms of energy 
loss to the phonon bath and discuss the proper modifica- 
tions of the master equation later in Sec. IVII 

We now apply secular approximation to Eq. (|16|) . This 
is justified due to the separation of the characteristic time 
scales Ej 1 <C r~ ut <C Tr 1 established in the previous sec- 
tion. When considering the evolution of the off-diagonal 
elements of the density matrix P + _(E p , t), we need to keep 
only terms oc P + _(E p , t) in the right-hand-side of the cor- 
responding master equation. The contribution of other 
elements of the density matrix to the evolution of the 
coherences P + _(E p ,t) is small as r out /Ej and T in /E r . 
Thus, we arrive at the equation governing the evolution 
of the coherences 

P + _{E p ,t) = -iE J P + _{E p ,t) 

i W +( E p> E k) + W-(E P , E k )] P + _(E P , t). 

k 

(17) 

The transition rates W + (E p , E k ) and W.(E p ,E k ) are 
given by the Fermi golden rule 

W+(E p ,E k ) = 2n\ (E P1 +\H T \o,E k ) \ 2 5(E p +LU+-E k ) 
W.(E p ,E k ) = 27r\ (E P ,-\H T \o,E k ) \ 2 5(E p +u_-E k ) 

(18) 

with uj± of Eq. (J2J). The matrix elements 
(E p , ±| Ht |o, E k ) can be calculated using the particle 
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conserving version of the Bogoliubov transformation 17 
and are (E p ,±\H T \a,E k ) = ^/2(t pk u p u k - t kp v p v k ) 7 
where Up , Vp HXC Bogoliubov coherence factors: 



„J = I(l + |),a»d„= = i(l-| 



(19) 



leading t< 



.18 



W + (E p ,E k ) = 2n\t pk 



2 (l + ^^fy(E p +oj + -E k ) 



W. (E p ,E k ) = 2n\t pk \ 2 (l + ^ 8{E P +U--E k ), 

(20) 

Now we may relate the tunneling matrix elements to the 
normal-state junction conductance 

9 T =8n 2 Y,Kk\ 2 S^p)S(Ck)- 

p,k 

Assuming that tunnel matrix elements t pk are weakly 
dependent on the energies £ k , £p, we can rewrite Eq. if^Hjl 
in terms of the dimensionless conductance: 



W.(E p ,E k ) = 



4~ 

g T S r Sb 

47T 



E p E k J 



(21) 

where 5 r r{,\ is mean level spacing in the reservoir (box) 
^r(6) = ( v r(b)Vr(b)) with v r (p) being single-spin electron 
density of states at the Fermi level in the reservoir (box). 

The system of equations for the diagonal part of the 
density matrix describes the evolution of the populations 
and follows from Eq. (|16|l . From now on we adopt the 
short-hand notation for the diagonal elements of the den- 
sity matrix (s\ p(t) \s) = P ss (E s ,t). In particular, we de- 
note the probability of the qubit to be in the state |+) 
or in the state |— ) and a quasiparticle to have energy E p 
as P ++ (E p ,t) or P (E p ,t), respectively; the probabil- 
ity P (E k ,t) corresponds to the state with a quasiparticle 
residing in the CPB and having energy E k . In these nota- 
tions the system of equations describing the dynamics of 
the populations for the states \+,E p ), \—,E p ) or |o,J5fc) 
can be written as 

P ++ (E p ,t)+J2w + (E p ,E k ) [P ++ (E p ,t)-P Q (E k ,t)] = 0, 
k 

(22a) 

P..(E p ,t)+J2W- (E p , E k ) [P . (E p ,t)-P o (E k ,t)]=0, 

k 

(22b) 

P (E k ,t)+J2 [W+(E P , E k ) + W-(E P , E k )} P Q (E k ,t)- 
p 

[W + (E p ,E k )P + JE p ,t) + W_(E p ,E k )P _(E p ,t)] = 0, 

p 

(22c) 



where we neglected the contribution of the coherences. 
This is justified as long as P + _(E p ,0) = in the initial 
moment of time, and the two parameters r ou t/-E 7 and 
Ti n /Ej, are small. The transition rates in Eqs. (|22|l are 
given by the Fermi golden rule, see Eqs. l|18fl . 

At the end, experimentally observable quantities can 
be obtained from Py by taking the proper trace over the 
quasiparticle degrees of freedom 



(23) 



where i, j = +, — . This completes the derivation of the 
master equations without quasiparticle relaxation, and 
we proceed to the solution of these equations. 



IV. EVOLUTION OF THE QUBIT 
COHERENCES. 

We now discuss the solution for the off-diagonal el- 
ements of the density matrix. We assume that ini- 
tially the qubit and quasiparticle are independent; the 
quasiparticle is in thermal equilibrium in the reservoir, 
and the qubit is prepared in a superposition state with 
a + _(0) ± 0: 



P + _(E P ,0) = Podd(E p )a + _(0) ? 0, 



(24) 



where p dd(E p ) is the equilibrium distribution function 
with an odd number of electrons in reservoir at temper- 
ature T < A, 



Podd(E p ) 



exp (-Ep/T) 



Z 



dd 



(25) 



The normalization factor Z Q dd here is 



Zodd = ^ exp 



E 



ttA T 

2j:VA exp v T 



The solution of Eq. (|17|) is straightforward. After trac- 
ing out quasiparticle degrees of freedom we obtain 

a + it) = a + S0)J2podd(E p )e^(-iE J t-lr ia (E p )t\,(26) 



where Ti n (E p ) is given by 

T in (E P ) = W+(E p ,E k ) + W.(E p ,E k )] . (27) 
k 

In the low-temperature limit T <C W- , w+ , the expression 
for <7 + _(t) can be simplified 



a + _{t) = a + _(Q)^p(-iE J t-^\ . 
Here the phase relaxation time T2 is given by 



8tt 



2A + w 



2A + u 



(28) 



(29) 
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The decay of qubit coherences is determined by the 
rate of the quasiparticle tunneling into the box as previ- 
ously discussed in Sec. |n] This result remains valid also 
in the presence of quasiparticle relaxation. On the con- 
trary, the evolution of the diagonal parts of the density 
matrix P++ and P— — , depends strongly on the relaxation 
of the quasiparticles. We will study this evolution with 
and without quasiparticle relaxation in the next sections. 



V. KINETICS OF THE QUBIT POPULATIONS 
WITHOUT QUASIPARTICLE RELAXATION 

The evolution of the diagonal elements of the density 
matrix is described by Eq. I|22l) . We will assume that 
initially the qubit is prepared in the state |+), and the 
quasiparticle resides in the reservoir. As explained in 
Sec. ITT1 tunneling out of the box (~r out ) is much faster 
than tunneling in (~r; n ) due to the differences in the 
volumes of the CPB and the reservoir. In fact, for a 
sufficiently small box, l/r ou t is the shortest time scale 
in the system of Eqs. Q22J1 . Therefore, we may neglect 
term d t P (Ek, t) in Eq. I|22cfl : i.e., the value of P Q (E k ,t) 
follows instantaneously the time variations of P ++ (E p ,t) 
and P__(E p ,t). This greatly simplifies the system of 
equations for the populations. The solution for P (E k , t) 
in this approximation is 



Pa{E k ,t) 



J2 p W + (E p ,E k )P ++ (E p ,t) 
J2 p [W + (E p ,E k )+W_(E p ,E k )] 

52 p W_(E p ,E k )P_(E p ,t) 
V [W + (E p ,E k ) + W_(E p ,E k )] 



(30) 



After substituting this expression back into Eqs. I|22a|) 
and (|22b|) . we obtain effective rate equations for the 
qubit in the presence of an unpaired electron in the su- 
perconducting parts 

P ++ (E p , t) + 7+ (E P )P ++ (E p ,t) - 7+ (E P )P . (E p +Ej ,t) = 0, 



P- - (E P , i) + 7- ( e p)P- - (E P ,t) - 7_ (E P )P ++ (E p - E 3 ,t) = 

(31) 

with j ± (E p ) having the form 

rFM _V W + (E p ,E k )W.(E p ,,E k ) 
7+1 P> £%J: p >>lW + (E p ,,,E k )+W_(E p ,,,E k )y 

(F] _ V W.(E p ,E k )W + (E p ,,E k ) 

{ P> P%Y.AW + {E p „,E k ) + W.{E pll ,E k )]- 



(32) 

This structure of the transition rates reflects the nature 
of the transitions involving an intermediate state |o, E k ). 
The normalization condition 



Y,[P ++ (E p ,t)+P.(E p ,t)} = l 



(33) 



is preserved under evolution. This can be checked di- 
rectly with the help of Eqs. I|31|l and the following rela- 
tion for the rates: 



7 + (E P )X(E p ) = Y."t- ( E v)X{E p - E, 



(34) 



(here X(E p ) is an arbitrary function of E p ). 

Let us discuss the solution of the Eqs. I|31|l . In the 
initial moment of time the qubit and quasiparticle are 
uncorrelated; the qubit is prepared in the excited state 
|+) and quasiparticle can be described by the equilibrium 
distribution function p dd{E p ): 

P ++ (E p ,0)= Podd (E p ) and P_(E p +E J7 0) = 0. (35) 

Upon solving Eqs. I|31|l . we find expressions for the pop- 
ulations of the qubit levels 



+ E 

p 

.(*) = £ 



Podd(E P h_ {Ep + Ej) 
1 _{E p + E J )+ 1+ {E p ) 

Podd(E P h + (E p ) 
1 _{E p + E J )+ 1+ {E p ) 

Podd(E p -Ej)j_ (E p ) 
1- (E p )+ 1+ {E p -E J ) 

Podd(Ep-E J )-f_(E p ) 
lAE P )+l + {E p -E J ) 



exp (-T(E p )t) , 



exp{-T(E p -E,)t) 



(36) 



where r(£" p ) is defined as 

T(E p )= 1+ (E p ) + 1 _(E p + E J ). 



(37) 



In the low-temperature limit, we calculate the sums in 
Eq. ijSHl) assuming A > u + > E r > T to find 



E 



l+(E p )p odd (E p ) 

nE P ) 



1 T 
ttE, 



(38) 



and 



a__(t) 



I T 
nE, 



T 
^E, 



■ exp 



t 

"f* 



I T 
ttE, 



1 — exp 



Tt 



The relaxation time Tj* is defined as 

9t S ' 



1 



4tt V 2A 



1 



CJ 4 



(39) 



(40) 



In deriving this expression we assumed E r / 'A <^ 1 and 
kept only the leading terms. 
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The solution for the populations in this case (no quasi- 
particle relaxation, r = oo) show that final qubit popu- 
lations are determined by the tunneling rates, which, in 
turn, depend on the superconducting DOS at different 
energies v(E p ) and v(E p + E 7 ). Since the states with 
higher DOS are more favorable, the quasiparticle can be 
found most of the time with energy close to A and rarely 
with energy A + E r . Therefore, at low temperatures 
T <C Ej, the qubit will mostly remain in the excited 
state |+). The probability to find the qubit in the ground 
state is proportional to ~ \JT jE 7 and thus is small, see 
Sec. ID] As soon as we include mechanisms of quasipar- 
ticle relaxation into consideration, the qubit populations 
will eventually reach equilibrium. In the next sections we 
investigate the equilibration of the qubit. 

VI. KINETICS OF THE QUBIT POPULATIONS 
WITH QUASIPARTICLE RELAXATION IN THE 
RESERVOIR 

A. Master equations with quasiparticle relaxation 

Here we consider a more realistic model by incorpo- 
rating the mechanisms of quasiparticle relaxation into 
the rate equations. Such mechanisms were studied in 
the context of non-equilibrium superconductivity^^. 
In aluminum, a typical superconductor used in charge 
qubits, the dominant mechanism of quasiparticle relax- 
ation is due to inelastic electron-phonon scattering. The 
relaxation time depends on the excess energy e of a quasi- 
particle^ 

J_ = 1 64^/A\ : W/ 2 

r(e) tq 105 \T C ) \A/ ' 1 j 

where e = E p — A <C A, and To is characteristic pa- 
rameter defining electron-phonon scattering rate at T — 
T c (here T c is superconducting transition temperature). 
For typical excess energies of the order of E, ~ 0.3 
K, the estimate for r yields quite long relaxation time 
r~10- 5 -10- 4 s. 

The procedure developed in Sec. IIIII allows us to in- 
clude the mechanisms of quasiparticle relaxation into the 
master equations. One can start by writing an equation 
of motion for the density matrix that includes the qubit, 
quasiparticle, and phonons, then expand the density ma- 
trix in the small coupling parameter - electron-phonon 
interaction as discussed in Sec. IIIII Finally, one should 
trace out phonon degrees of freedom and obtain master 
equations for the qubit with quasiparticle relaxation. We 
will skip the cumbersome derivation and present only the 
results here. In the relaxation time approximation the 
collision integral has the form 

l±=~(P ±± (E p ,t)-P ±± (E p ,t)), (42) 

where t — t(e^Ej). The probability P ±± (E p ,t) is pro- 
portional to the equilibrium distribution function of a 



quasiparticle p dd(E p ) and the proper qubit population 
a ±± (t): 

P ±± (E p ,t) = Podd (E p )J2P ±± (E p ,t). (43) 
p 

The form of P ±± (E p ,t) is dictated by the fact that 
phonons equilibrate the quasiparticle only, without af- 
fecting directly the qubit state* 1920 . The collision in- 
tegral Eq. H42fl replaces zero in the right-hand sides of 
Eqs. l(22a|) and $Zty . However, Eq. (|23c|) for P {E k ,t) 
remains unchanged due to the short dwelling time of a 
quasiparticle in the box (we assume that r(e ~ lo + ) 
r~ ut , but set no constraints on rFi n ). Then, the system 
of Eqs. (|31|l for populations can be written as 

P ++ (E p , t) + 7+ (E P )P ++ (E p , t) - 7+ (E p )P..(E p +E Jt t) = 
-±(P ++ (E p ,t)-P ++ (E p ,t)), 

P_ _ (E p , t) + 1 _ (E p )R . (E p ,t)- 1 _(E p )P ++ (E p - Ej , t) = 
-^(P.(E p ,t)-P__(E p ,t)) (44) 

with the effective transition rates 'y ± (E p ) defined in 
Eq. (|32|l . The obtained system of integro-diffcrcntial 
equations l|4^|) for P ±± (E p , t) describes the effect of quasi- 
particle relaxation on the dynamics of the qubit. 

We solve Eqs. I|44|) first in the simple case of a short 
relaxation time (r rr ). Under these assumptions, 
we can seek the solution in the form 

P ±± (E Pl t) =p odd (E p ) a ±± (t), (45) 

with <J ±± (i) defined in Eq. so that P ±± (E p ,t) = 

P ±± (E p , t). Using this ansatz and performing the appro- 
priate summation, Eqs. I|44|) reduce to the Bloch-Redfield 
equations 

a ++ (*) + ^2l+( E p)Podd(E p )a ++ (t) - 
p 

p 

*--(*) + ^7-(£p)podd(£ P )M*)- 
p 

- ^7_(£ P Kdd(£p -£>++(*) = o. 
p 

(46) 

Utilizing the property of the rates Eq. (|34|) , one can sim- 
plify the equations above, 

*++(*) + <7 + >°- ++ (*) = <7->o---(t), 
M*) + <7-><Mt) = ( 7+ )(7 ++ (t), (47) 
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where thermal-averaged transition rates ^7 + ) and (7 ) are of algebraic equations 

<7+) = ^2l + (Ep)Podd(Ep), 
v 

<7_>=I>-(£ P )Podd(£y. (48) 



One can also check that due to relation l)34|l rates (7 ± ) 
comply with the detailed balance requirement 



<7 + ) 



/ \ _ eXP \ T 

(7_) V T 



For the initial conditions: er ++ (0) = 1, <7__(0) = 0, the 
solution for populations is 



«■++(*) = 



a..(t) = 



-EJT 



,-((7 + ) + (7_))* 



l + e- E .i/ T l + e- E j/ T 



1 



1 - e 



-((7 + )+(7_))* 



(49) 



In the low-temperature limit (T -C i?, < w + < A) we 
find a simple form for the effective rates (7 ± ) in the lead- 
ing order in T/A and E T /A: 



<7-> 



4tt V 2A 



4tt V 2A + 



/ r 

7rEj 

ry 

ttE, 



■ exp 



T 



(50) 



Factors of \JT jitE } in the rates can be interpreted as 
the probability of flipping the qubit (|+) — * |— )), which 
is mainly determined by the ratio of the DOS of quasi- 
particles at energies v{E p ) and v{E p + E } ), respectively 
(see the discussion in Sec. [HJ|. We would like to point 
out here that energy relaxation rate (7 + ) + w_) is smaller 

than the phase relaxation rate by a factor of ^/T/ttEj. 



B. General solution for the qubit populations in 
the relaxation time approximation 

In this section we find the solution of Eqs. Ij44(l at an 
arbitrary value of T[ n T. In order to find the solution for 
the qubit populations we will use Laplace transform 



/>oo 

P{E p ,s) = / dtP(E p ,t)( 
Jo 



(51) 



and reduce the system of differential equations (|44|l sup- 
plied with the initial conditions Eq. (|35[1 to the system 



sP ++ (E p ,s) - Podd (E p )+ 7+ (E p )P ++ (E pi s)- 1+ (E p )R.(E p , S ) 

= ~{P ++ (E p ,s)-P ++ (E p ,s)), 
sP..(E p , s) + j_(E p )P_(E p , s)-7_ (E p )P ++ (E p , s) = 

= ~(?..(E p ,s)-P..(E p ,s)y (52) 



Here tilde denotes the shift by E f of the energy argu- 
ment in a function, e.g. P..(E p ,s) = P..(E p + E ; , s). 
The system of algebraic equations lf52*f can be solved for 
P ++ (E p ,s) and P {E p ,s). Then, by summing these ex- 
pressions over the momenta p and utilizing Eq. I|43l) we 
obtain a closed system of equations for qubit populations 
cr ±± (s): 

, \ ., . A(s) . . B(s) . , 
a ++ (s) = A(s) + -^a ++ (s) + -^o--(s), 

o-..(a) = C{s) + ^la ++ (s) + ^M CT __( S ), (53) 

where the coefficients A(s), B(s), C(s) and D(s) are 
given below 



s + 1/t 



s + 1/t ' 



s + 1/t 
Function -Z(s) is 



s + 1/r 



7M _ 7 + (-gp)Podd(^p) 

Z(sj ~^ s + i/T + r(^) 



(55) 



with jAE p ) and r(£7 p ) being defined in Eqs. J2TJ>, 
and (|37[) . respectively. From now on we take thermo- 
dynamic limit and replace the sum by the integral in 
Eq. (|55|) . (Thermodynamic limit is appropriate here, 
since S r <C T.) 



Z { s) = ^r d E pV{ E p f^}^ 
y ' S r J A p K v> s + l/r + TiE, 



(56) 



The solution of Eqs. i|53|) yields the following results 
for cr ±± (s): 



(ts + 1)Z(s) 



a ++ (s) = - 1 



+ Z(s)(l + e- E J / T ) 



Ms) = l_ (rs + l)Z(s) 



S TS 



+ Z(s){l + e- E j/ T )' 



(57) 



Equations l|57|) allow us to analyze the dynamics of the 
qubit populations for arbitrary Fi n r. Let us point out 
that Eqs. H57fl satisfy normalization condition er (s) + 
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<7__(s) = 1/s. To find the evolution of the populations, it 
is sufficient to evaluate a ++ (t). 

The inverse Laplace transform is given by 

= dsa ++ (s)e st , (58) 

L ~ % Jri-ioo 

where rj is chosen in such way that cr ++ (s) is analytic 
at Re[s] > rj. The integral (|58|l can be calculated using 
complex variable calculus by closing the contour of inte- 
gration as shown in Fig. |3 and analyzing the enclosed 
points of nonanalytic behavior of u ++ (s). In general, 
the singularities of <r ++ (s) consist of three poles and a 
cut. The latter is due to the singularities of the func- 
tion Z(s) causing er ++ (s) to be nonanalytic along the cut 
s G (s min , s max ), where 

Smin = max [T(E P ) ] , 

r 

Smax = mm[T(E p )] . 

T 

The schematic plot of T(E P ) is shown in Fig.^J 

In addition to the cut, a ++ (s) has 3 poles. The first one 
is at si = 0; two more poles, s 2 and S3, are the solutions 
of the following equation in the region of analyticity of 
the function Z(s): 

ts + Z{s) (l + e~ E J /T ) = 0. (59) 

The preceding discussion of the analytic properties of 
<7 ++ (s) is general for any ratio of the relaxation time r 
and quasiparticle escape rate r m . However, the location 
of the singularities and their contribution to the inte- 
gral depends on r m T. Below we briefly present re- 
sults for two cases of interest: fast (r in r -C 1) and slow 
(Ti n r 3> 1) quasiparticle relaxation. The detailed analy- 
sis of the singularities of er ++ (s) is given in the Appendix. 

In the fast relaxation regime (r; n T <C 1), the contribu- 
tions from the cut and the residue at S3 are small (pro- 
portional to ri n r, see the Appendix) and thus can be 
neglected. Then, relevant poles of a ++ (s) in this limit 
are 

si = 0, s 2 = - < 7+ ) - ( 7 _) . (60) 

The integration of Eq. i|58[) yields, up to corrections van- 
ishing in the limit T^t <C 1, Eqs. I|49|) for the popula- 
tions. 

In the slow relaxation case, r m r 3> 1, the main contri- 
bution to the integral l|58|) comes from the cut, and poles 
si = and s 2 . The latter may be found by iterative 
solution of Eq. 

s 2 = -M (l + e- V) , 

where function Z(jS) is defined in Eq. H5(jf) . and in the 
case of slow relaxation can be approximated as 

1 f°° t (E ) 

Z{Q) ^5~J A dE pHE P )f^yPo M (E p ). (61) 



Im(s) 





= — e e-< 


- 






mn S nmx S J S i S l 




Re(s) 



FIG. 3: (color online). Contour of integration (red line) cho- 
sen to calculate inverse Laplace transform Eq. 1581 . Points of 
nonanalytic behavior of <j ++ (s) are shown. (Poles at si, S2, 
and s 3 , and a cut s £ (s min , s max )). 

The residue at S3 gives a smaller by factor l/Ti n T <C 1 
contribution, see the Appendix. 

Taking integral in Eq. i|58|) along the contour enclosing 
the cut shown in Fig. |3 and accounting for poles Si and 
s 2 , we find 

e- E .i/ T ( 1 \ 
xexp(-Z(0)(l + e- B >/ T )^) + 

S~ J a LLE ) exp(~T(E p )t) . 

(62) 

Here we neglected the corrections to Eq. I|62|l of the order 
i/Lm T - The obtained expression forcr (£) describes the 
kinetics of the qubit populations in the slow relaxation 
regime. Note that the solution Eq. I|62|l satisfies initial 
conditions a ++ (0) — 1 and is consistent with previous 
results. Indeed, in the limit t— >oo, the exponent in the 
second term goes to zero and we recover Eq. (l3*fj|) (To 
show this one should use Eq. l|34fl L 

Equation (|62|l becomes physically transparent in the 
low-temperature limit. Using the approximation (|38l) for 
Z(0) and for the integral in the last term of Eq. (|62|) we 
find 
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where relaxation times Ti and T* are 



10 



in 



10 1 



The same quasiparticle density 




and 4 = ^ 



4tt V 2A 



1 + 1(64) 



Obtained results describe the relaxation of the qubit 
populations in the slow relaxation limit (r in r 3> 1) dis- 
cussed qualitatively in Sec.|nJ According to Eq. I|tj3|) . in 
this case the process of equilibration of qubit populations 
occurs in two stages. The first stage (t ~ ) corresponds 
to a quasistationary state formation with qubit popula- 
tions much larger than equilibrium ones. For the typi- 
cal experimental temperatures T ~ 20 mK, the excited 
state population of the qubit is about 85%. The equi- 
librium populations are established in the second stage, 
on the time scale of T\. The relaxation time Ti sets an 
important experimental constraint on the frequency of 
repetition of qubit experiments. 

We estimate now energy and phase relaxation times 
for the realistic experimental parameters 1 ^: A i=s 2K , 
E c « 0.5K, Ej w 0.3if, g T ~ 1, and T w 20mK. For the 
volume of the reservoir 10 -19 — 10 _17 m 3 , Eq. (|29[) yields 
the phase relaxation time T2 ~ 10~ 5 — 10 _3 s. Energy 
relaxation depends on the relation between Ti n and 1/r. 
Taking t~ 10~ 4 s and F" 1 ~ 10 _5 s, which corresponds to 
the lower end of the volume range, energy relaxation is 
described by Eq. f2J| with Ti*~10~ 5 s and Ti~10~ 3 s. 

We considered so far the effect of a single quasiparti- 
cle on the qubit kinetics. It is possible to generalize our 
results onto the case of many quasiparticles N qp residing 
in the system. In this case Ti becomes shorter since the 
quasiparticle tunneling rate r; n , see Eq. ©, should be 
multiplied by the number of quasiparticles in the super- 
conducting reservoir 



1 

T~ 2 



Tin • N, 



qp 



(65) 



(Here v is the normal density of states per unit volume 
and n qp is the density of quasiparticles in the reservoir.) 
Note that the volume of the reservoir does not enter in 
Eq. I|tj5|) ■ A finite density of quasiparticles in the reservoir 
affects also the process of energy relaxation of the qubit. 
The most clear example corresponds to the limit r; n r <C 1 
in which quasiparticle relaxation in the reservoir occurs 
fast compared to the time needed for the quasiparticle to 
reenter the CPB. In this limit T 1 /T 2 ~ y/Ej/T. 

The comparison of the theoretical prediction for T2 
and Ti with experimental data is complicated by the un- 
known value of iV q p in a qubit. The quasiparticle density 
in a system with a massive lead is known to deviate from 
the equilibrium value in a number of experiments^^. 
The estimate of n qp can be obtained from the kinetics 
of "quasiparticle poisoning" studied in the recent exper- 
iments 2324 . The observed rate of quasiparticle entering 
the CPB was 10 5 - 10 4 Hz. Assuming that the bot- 
tleneck for the quasiparticles was tunneling through the 
junction (rather than the diffusion in the lead), we es- 
timate the density of quasiparticles in the lead to be 



in a qubit with the reservoir volume 10 



-19 



would result in N, 



qp 



1 - 100. 



VII. CONCLUSIONS 

We studied the kinetics of a superconducting charge 
qubit in the presence of an unpaired electron. The pres- 
ence of a quasiparticle in the system leads to the decay 
of quantum oscillations. We obtained master equations 
for the coherences and populations of the qubit, which 
take into account energy exchange between the quasi- 
particle and the qubit, and include the mechanisms of 
quasiparticle relaxation due to electron-phonon interac- 
tion. Finally, we found decay exponents governing the 
dynamics of the qubit for different cases: fast and slow 
quasiparticle relaxation in the reservoir. 

We have shown that phase relaxation is determined 
by the quasiparticle tunneling rate to the box r in ~ 
g T 5 r /4ir. Kinetics of the qubit populations depends on 
the ratio of the quasiparticle relaxation time r and es- 
cape time rr . In this paper, we considered two limits 
- fast (rri n <C 1) and slow (rr in 3> 1) quasiparticle re- 
laxation. In the latter case, decay of qubit populations 
occurs in two stages. In the first stage at t ~ Fr 1 a 
quasistationary regime is established with large nonequi- 
librium excited state population. The second stage de- 
scribes the attainment of the equilibrium populations and 
occurs on the time scale of t^/ttE ~ /T. In the fast relax- 
ation case, equilibrium qubit populations are established 
att^T^^EjT. 
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APPENDIX A: ANALYSIS OF THE 
ANALYTICAL STRUCTURE OF a ++ (s). 

In this appendix we study the analytic properties of 
<7 ++ (s) in order to calculate the inverse Laplace trans- 
form (|58|l . In general, the nonanalytic behavior of er ++ (s) 
is determined by three poles, one of them is at s = 0, and 
a cut as shown in Fig. El The locations of two other poles 
and of the cut, and also contributions of all the mentioned 
singularities in a ++ (s) to the integral (|58p. depend on the 
value of r; n T. 

In the fast relaxation regime (Ti n T <C 1), in the vicinity 
of the s = pole, we find 



~ E .,/ T l 



-BIT 



(Al) 
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Two other poles S2 and s 3 are the solutions of Eq. ljl)9"|) 
with s 2 being the solution at small s ~ Ti n and s 3 at 
large s ~ 1/r: 



s 2 = - «7 + ) + <7-» , 

S3 = ---min[r(^)] + %^Ci. 

T ATT 



(A2) 



Here r(£ ; p ) is defined in Eq. l|37|) and C\ is a distance 
from the beginning of the cut in units of g T S r /4ir, see 
Fig. |21 In the vicinity of the second pole, cr ++ (s) is given 
by 



l + e- E j' T 



Z(0) 



32 



(A3) 



with Z(0) defined in Eq. ijfiTJl. The residue of c ++ (s) at 
s 3 is proportional to (Tr in ) 2 . Consequently, the contri- 
bution to the integral (|58|l from the pole s 3 is small. 

In addition to the poles discussed above, nonanalyt- 
icity of cr ++ (s) comes from the singularities of Z(s). 
The function Z(s) is nonanalytic along the cut s £ 
(s m in,s max ), where 



Smin = ---max[T(E p )] 



Smax = min[r(£L)] , 

T 



with T(Ep) being defined in Eq. I|37|l . The proper contri- 
bution to Eq. i|58|) can be calculated by integrating along 
the contour enclosing the cut 



-1 
2ttI 



dse st (cr ++ (s+ie) — <J ++ (s— iej) 



(A4) 



The discontinuity of the imaginary part of cr ++ (s) at the 
cut is 

<J ++ (s+ie) — cr ++ (s — ie) = 

-2iT(TS+l)ImZ(s+ie) 
~ [Ts+ReZ(s)(l + e- E j/ T )] 2 + [lmZ(s + te)(l + e" E j/ T )] 2 ' 

(A5) 



In the limit Ti n T -C 1 we find 



<T ++ (s + ie)-<T ++ (s — ie) 



-2iT(TT ln )lmZ(s + ie) 



[-l+ReZ(s)] 2 + [lmZ(s+ie)} 2 ' 

(A6) 



which yields a negligible contribution to Eq. I)58|l from 
the cut, I cut oc TinT <C 1. Finally, after summing up two 
relevant contributions, one obtains the result for a ++ (t) 
given in Eq. ig^J. 

In the opposite limit of slow relaxation (r^r 3> 1), the 
first pole s% = is the same as in the previous case with 



m 

r 







r( E ) 






r(0) 








\ r n,i„ 










£ 



FIG. 4: (color online) Main panel: dependence of T(e), de- 
fined in Eq. jjffi ■ on quasiparticle energy e = E p — A. Inset: 
dependence of T(e) in the vicinity of £ min . The maximum 
value of the rate is r max = g T 8 r /in (dashed line). In the limit 
E j <C uj + <C A, a simple analytic expression for the minimum 

can be found: e min = _E 7 /25, and rmi °~ r (°) 



r(o) 



C— , where 



r(0) = 4r t, v/ ^4-, and numerical constant C « 0.1. 



the expression for <r ++ (s) given by Eq. (jJTJ. The other 
two poles S2 and s 3 are found from Eq. (|59")l assuming 
r in r » 1: 



s 2 



■S3 



-min[r(£L) 



Co 



T 2 T- 



(A7) 



where Z(0) is defined in Eq. (|61|l . and C2 is a positive 
constant of the order of unity. In the vicinity of the 
second pole <r ++ (s) is given by 



C++0) = 



1 



l + e~ E j/ T 



Z(0) 



1 



S2 



(A8) 



The third pole s 3 lies in close proximity to the begin- 
ning of the cut s max . In order to find the position of the 
pole we expand the denominator of Z{s) in the neighbor- 
hood of the minimum of T(E p ), see Fig.0] 



Z(S)- 



dE„ 



v{E P h + (E p )p odA {E p ) 



A S r a + l/T+r aSn +r"(E Pjnin )(E p -E Pm Js 



(A9) 



and solve Eq. (|59|l to obtain s 3 of Eq. (|A7fl . Here we 
used the small-?/ asymptote, Z(y) oc 1/^/y, where y is 
the distance from s max in units of T^; i.e., we made a 
substitution 



= -i/r-r„ 



Tin?/ 



in Eq. (|A9|> and evaluated the integral at y <C 1. Then, 
in the neighborhood of s 3 the expression for tr ++ (s) can 
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be written as 

cr ++ (s) « 



r in rZ(j/ 



r in r + Z'(yo)(l + e-^/ T )s- S3 : 



(A10) 



where ?/o is the corresponding solution of Eq. (|59|l and 
is small: yo ~ (Fi n r)~ 2 <C 1. Since derivative Z'(yo) cx 

j/ 3 ^ 2 , the residue at S3 is suppressed as (TjnT) -1 and 
thus can be neglected. 

The contribution from the cut in the limit T[ n T 3> 1 
can be evaluated from Eq. (| A5|) . The discontinuity of the 
function cr ++ (s) is 



c ++ (s +ie) - <7 ++ (s — ie) 



2i 



lmZ(s+ie). (All) 



Hence, the contribution to Eq. (|58|l from the cut is 



lent — — 



ds—S(s+T(E p )) 



00 dE 

-j^v(E p yf + (E p )p odd (E p ) 

°° dE p ^ + (E p ) Podd (E p ) p/px., / a10 n 

-^^(^p)— - — \ expl-T(E p )t) . (A12) 



Finally, combining proper terms one finds the inverse 
Laplace transform of Eq. (|62() . 
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